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Abstract 

In this paper we investigate some entanglement properties for the Hy- 
drogen molecule considered as a two interacting spin J (qubit) model. 
The entanglement related to the H2 molecule is evaluated both using the 
von Neumann entropy and the Concurrence and it is compared with the 
corresponding quantities for the two interacting spin system. Many as- 
pects of these functions are examinated employing in part analytical and, 
essentially, numerical techniques. We have compared analogous results 
obtained by Huang and Kais a few years ago. In this respect, some pos- 
sible controversial situations are presented and discussed. 

1 Introduction and the model 

Entanglement is a physical observable measured by the von Neumann entropy 
or, alternatively, by the Concurrence of the system under consideration. 

The concept of entanglement gives a physical meaning to the electron cor- 
relation energy in structures of interacting electrons. The electron correlation 
is not directly observable, since it is defined as the difference between the ex- 
act ground state energy of the many electrons Schrodinger equation and the 
Hartree-Fock energy. 

In this paper we discuss the Hamiltonian which describes the Hydrogen 
molecule regarded as a two interacting spin | (qubit) model. 

In [1] it was argued that the entanglement (a quantum observable) can be 
used in analyzing the so-called correlation energy which is not directly observ- 
able. From our point of view, the Hydrogen molecule is dealt with a bipartite 
system governed by the Hamiltonian 

Uh 2 = + 9)^1 ® cti - ^-(1 - g)o- 2 ® o- 2 - B(a 3 ® <t 3 + a <g> cr 3 ), (1) 
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where <jj stand for the Pauli matrices (<tq — I). Actually, this model was con- 
sidered in [I] in order to illustrate their method. However, here we will make 
some interpretative changes. Indeed, from our point of view, the states of an 
isolated atom are strongly reduced to a system with two energy levels related to 
the intensity of the magnetic field B. Relatively to this scale, the exchange in- 
teraction constant J is usually smaller than B, in order to represent the residual 
interatomic interactions. From the point of view of quantum chemistry, one may 
interpret the discrete spectrum as provided by the Hartree-Fock calculations, 
while the interaction coupling J models the residual multielectronic effects, not 
taken into account by the mean field approximation. 

For simplicity we limit ourselves to the ferromagnetic phase with J > 0. 
The parameter g, such that < g < 1, describes the degree of anisotropy 
corresponding for g — to the completely isotropic XY spin model. Conversely, 
g=l provides the anisotropic XY spin model, the so-called Ising model. 

We notice that when the atoms are far apart, their interaction is quite weak. 
This corresponds to a vanishing value of J. In this situation the state of the 
system is completely factorized in the product state of the ground states of the 
indipendent spins. The corresponding total energy, in unit of B, is just the sum 
of the two fundamental levels, £ a = —2, which we may consider as the Hartree- 
Fock approximated fundamental level in molecular structure calculations. 

When J 7^ 0, the fundamental energy eigenvalue is £= —y^4 + g 2 X 2 in Re- 
gion I defined by < A < . 2 , otherwise £ = — A (A means the coupling 

\/l-9 2 

constant) in Region II, which is the complement of I which respect to pos- 
itive real axis. The corresponding (non normalized) eigenstates are — 

{ ^ S gA + ~~' 0' Q> l} ana - l^n) = {o, 1,1,0 j, respectively. In both cases the 
state is entangled. 

Since we are dealing with pure states, the von Neumann entropy [2] 

SvN = -Tr[pilog 2 pij (2) 

is chosen to be a measurement of the entanglement, where p\ is the 1-particle 
reduced density matrix. However, for general mixed states other entanglement 
estimators (for instance, the Concurrence [3]) have to be used. In the considered 
case, one has 



S V N.\ — — 



. 9 2 vVA 2 +41og(l 



1 



vVA 2 + 4 

■ (d 2 (Vg 2 * 2 + 4 + 4) A 2 + 8 (V.g 2 A 2 + 4 + 2) ) log (i 



A 2 



2 VPA7T4 



1 



(3) 



vVA 2 +4 [g 2 \ 2 + 2 v VA 2 +4 + 4J log(4) 

S vN ,n = 1. (4) 
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Scrutinizing Eq. and Eq. (g]) it emerges that the entropy is an increasing 
function of the coupling constant A in Region I, but the state is maximally 
entangled in Region II independently from the anisotropy parameter g. One 
sees that, as it arises graphycally, for g = 1 the entanglement is a monotonic 
increasing function of the interaction coupling A. Moreover for weak (< 1) 
coupling values it is always less than the 30%. Of course, for large coupling- 
constants the entropy approaches 1, meaning that all levels are equiprobably 
visited by the considered spin. 

Limiting all further considerations to the case of weak interaction, we observe 
that at the boundary point Ab = . 2 a discontinuity occurs, signaling a 

crossing of the lowest eigenvalues and, in a more general context, a quantum 
phase transition [5]. 

As it was pointed out in [B], for quantifying the entanglement we can resort 
to the reduced density matrix. Furthermore, in [7], Wootters has shown that 
for a pair of binary qubits one can use the concept of Concurrence C to measure 
the entanglement. 

The Concurrence reads 

C{p) = max(0, V\ — v-x — v% — v±), (5) 
where the iVs are the eigenvalues of the Hermitian matrix 

R = (VppVp) 2 ) 

where p — (a y ® a y )p*(a y ® o y ), p* being the complex conjugate of p taken in 
the standard basis [7]. 

Some interesting results on the simple model |T]) of the Hydrogen molecule 
can be achieved by realizing a comparative study of the von Neumann entropy 
and the Concurrence. 

To this aim, we compute the Concurrence C\ and Cn, i- e. 

where I and II refer to Regions I and II, where < A < j^i, and £ = — A, 
respectively. 

In Figure [1] a comparison between the Concurrence and the von Neumann 
entropy for two spins system as a function of the coupling A for g = 1 is pre- 
sented. 

Sec. [5] contains a comparison between the entanglement and the correlation 
energy. In Sec. [3]the Configuration Interaction method is introduced to compare 
entanglement and correlation energy. In Sec. [4] some differences between the 
Configuration Interaction approach and the two spin Ising model are presented. 
Finally, our main results are summarized in Sec. [5) 
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2 A comparison between the entanglement and 
the correlation energy 

Now we look for a comparison between the entanglement with the energy cor- 
relation, which as we have already recalled, it is understood as the difference of 
the fundamental energy level compared with respect to the corresponding value 
at vanishing coupling constant A. 

For g — 1 and in unities of B it is given by 



£corr = \£o\ ' 2 = V4 + A 2 - 2. (7) 

We observe that the entanglement measure is always bounded, while £ cor r is 
a divergent function of X. So it does not make much sense to look for simple 
relations valid on the entire X-axes. Consequently, limiting ourselves to weak 
couplings, for < A < 1, we minimize the mean squared deviation 
~i 

AS^dX, with AS a = £ cor r — a. S v n ■ (8) 



a 


-I 







Thus the minimizing parameter a m i n will be given by 

Jo £corrS v N dX 



Jo S% N dX 



0.691217. (9) 



A formula analogous to © can be obtained by using the Concurrence as 
a measure of entanglement. In this case, by minimizing the mean squared 
deviation we have 

= f AC^dX, With AC a > = £ C orr — Ot C. (10) 

Jo 

Now, in order to estimate the relative deviation of S v n with respect to £ CO rr, 
let us report |AS' Qmira |/iS' t ,^r and \AS amin /£ corr \ as functions of A at the optimal 
value amin- The graphs of these functions are shown in Figure [2] 
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AS ml „/S vN | |AS mln /E corr | 




A 
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Figure 2: The relative quadratic deviation between the von Neumann entropy 
and the correlation energy with respect to the former and the latter, respectively, 
at the optimal value a m in as a function of the coupling constant A for g = 1. 

In Figure [3j the relative quadratic deviation between the Concurrence and 
the correlation energy with respect to the former and the latter, at the optimal 
values a' min , is represented. 



AC ml „|/C |AC min /E corr | 
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Figure 3: The relative quadratic deviation between the Concurrence and the 
correlation energy with respect to the former and the latter, respectively, at the 
optimal value a' min as a function of the coupling constant A for g = 1. 

Remark 1 

From these graphs, one can argue that the agreement between the two func- 
tions S v n and Ecorr is only qualitatively good, in fact, for very small A, it is not 
good at all. However, in an intermediate range of values, i. e., 0.6 < A < 1 the 
two functions are almost proportional within the 10%. Analogously, the same 
is true between energy and Concurrence. Even, the agreement becomes worst 
comparing the relative deviation of the Concurrence with respect to the corre- 
lation energy, since the range in which the relative deviations become smaller 
than 10% are narrower. Then, the question is whether the above results are 
i) sufficient to justify the conjecture advanced in pQ, i.e., entanglement can be 
considered as an estimation of correlation energy; ii) if such a relation has a 
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more concrete physical meaning, in particular whether the minimizing parame- 
ter a m in and the vanishing point of AS amin does possess any physical meaning 
(or a' min and the vanishing point of AC a > . ). Notice that in the case of the 
comparison for the Concurrence simpler analytical expressions appear. For in- 
stance one finds = (o.383249 ^ +4 - \A 2 + 4 + 2) . 

Remark 2 

We note that in an interval of values around a m in , the deviation function 
(J8j possesses a minimum in the interval of interest < A < 1, otherwise the 
minimum is achieved at larger value of A, or the function is monotonically 
increasing (see Figure 0J. 



AS a d AS a / dA 




Figure 4: The deviation AS a and its derivative with respect to A are computed 
for values of — 1.29(red) < a < — 0.091(violet), for steps of 0.06. The curve 
drawn thicker corresponds to a m in 



This behavior suggests to consider the function AS Qmjn as a sort of "free en- 
ergy" , where a m i n mimics the " temperature" specific of the system. If, for some 
reason, we allow A to change, then we expect that spontaneously the interaction 
coupling adjusts itself to the minimum of AS amin . Similar considerations can 
be made looking at the graphs drawn for the function AC a > and its derivative 
with respect to A (see Figure [5]). 

The function AS amin or, alternatively, the minimum of AC a > . can be ob- 
tained algebraically. Such a minimum is at the value of the coupling constant 

A m"in ~ °- 485 and A Ln ~ °- 371 ' respectively. 

The authors in pQ studied numerically the von Neumann entropy and the 
correlation function for a Hydrogen molecule, using an old result by Herring and 
Flicker \$. , going back to an oldest idea by Heitler and London [9] , which con- 
sists in substituting the molecular binding with a position dependent exchange 
coupling: 

J(r) « 1.641 ri e~ 2r Ry, (11) 

where r is given in Bohr radius, see Figure[6l The maximum value taken by this 
function is at the point r max — 1.25. Assuming B = 0.5 Ry, i.e. \ of the funda- 
mental level of the Hydrogen atom, the maximum value X' max = J(r max )/B ~ 
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Figure 5: The deviation AC a > and its derivative with respect to A are computed 
for values of — 0.98(red) < a' < 0.22(violet), for steps of 0.06. The curve drawn 
thicker corresponds to a' min 

J(r) (Ry) 




Figure 6: The effective interaction Hydrogen-Hydrogen atom 



0.470628 < A m "^, i.e. the value of the effective interaction value is less than 
the minimum for the deviation function AS amin . Then, the equilibrium bal- 
ance between entanglement (as von Neumann entropy) and correlation energy 
predicts a length of the molecule equal to r max (see the first panel of Figure [7]) . 
On the other hand, if we consider the energy gap 2B = 3/4 Ry, i.e. the energy 
step to the first excited state, one obtains the new value X'' x ~ 0.628, which 
goes beyond \ m im even if it is always less than 1. Now, the deviation function 
AS amill has two minima as seen in the second panel of Figure [71 one of which 
is at r"_ w 0.76 , the other one being at r'l « 1.91. 

These results should be compared with the experimental equilibrium length 
of the Hydrogen molecule, which is r exp ss 2.0. 

We point out that although the spin-model described by the Hamiltonian 
([T]) is characterized by features which are essentially rough, however we are 
induced to answer positively to the quest for a physical meaning of the deviation 
function A5 amin . Indeed, the results elucidated in Figure [7] encourage, on one 
part, improvement of the computation of r in order to make more accurate the 
comparison with the experimental value r exp . 
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AS mln (r; B= .5 Ry) AS^fr; B= .375 Ry) 




Figure 7: The von Neumann entropy for the 2-spin model for B = .5 Ry (left 
panel) and for B = .375 Ry (right panel) and the position depending interaction 
given by (fTTj) . 

The first question to answer is whether this draft works also for the Concur- 
rence. A statement about it is not obvious, since the von Neumann entropy is 
a nonlinear function of the Concurrence in the 2-qubits case. 

However, from Figure [8] one can see that the minimized deviation of the 
Concurrence takes one minimum for relatively large intensity of the magnetic 
field ( say B > 0.6 Ry), while for weak fields two minima appear, corresponding 
to the situation depicted nearby. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 8: Two contour plots of the minimized deviation of the Concurrence as a 
function of the magnetic field B (Ry) and of the internuclear distance r, as given 
by (fTTj) . The range of values divided by the contour lines is [—0.038, 0,04] for 
the left panel and [—0.03705, —0, 03000] for the right one that approximatively 
corresponding to the black area in the left panel. 

In correspondence of the same values considered above, for B = 0.5 Ry 
the function AC a > . (r) has two minima at r = 0.79 and r = 1.88, while for 
B = 0.375 Ry they are located at r = 0.60 and r = 2.25. So one sees that the 
resulting equilibrium configurations are not much very close to the experimental 
one. The equilibrium configuration more closest to the experimental one is the 
minimum occurring at r = 1.88 [B = ^ Ry) for the function AC a > (r). 
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One sees that one of the resulting equilibrium configurations is only roughly 
close to the experimental one. 

In other words, to conclude monitoring numerically B the equilibrium config- 
uration more closest to experimental one in the minimum occurring at r = 1.88 
for B = \ Ry and at r = 2.25 for B = 0.375 Ry for the function AC a > . (r). 



3 A quantum chemical framework to compare 
entanglement and correlation energy 

In this Section we represent the results produced in [T] , where the electron entan- 
glement in the Hydrogen molecule, calculated by the von Neumann entropy of 
the reduced density matrix p\ , is obtained starting by the excitation coefficients 
of the wave function expanded by a configuration interaction method: 



' ' " D log 2 pf ISD ^ 
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(12) 



where c\ is the coefficient for a single excitation, and ci^ is the double 
excitation (in Appendix A of [TU] more details are shown). 

In this framework, entanglement (S) and correlation energy (£ corr ), as func- 
tions of nucleus - nucleus separation are those in Figure [9] 





0.8 


to 




c 


0.6 


E 




angh 


0.4 


c 

UJ 






0.2 








H 2 






m ■8(p ] ) 1, 










y / 



















2 3 

R(A) 



0.05 
0.04 
0.03 



0.02 r 

UJ 



0.01 




Figure 9: Comparison between the entanglement, calculated by the von Neu- 
mann entropy of the reduced density matrix, and the electron correlation energy 
in the Hydrogen molecule. 



By the results given by this model, we want to discuss and to suggest some 
answers to the questions i) and it) presented in Remark 1. Even if, in order 



9 



to represent correlation energy and entanglement, we use two different scales, 
in Figure [9] we can sec that entanglement has a small value in the united atom 
limit after it is growing for small distances till it arrives at a maximum value 
then it decrease till it assumes zero value at the separated atom limit and it is 
exactly the progress of the correlation curve. 

In order to compare the entropy S with the electron correlation energy £ cor r , 
we rescale S with the parameter a m i n calculated with some procedure illustrated 
in Eq. ([5]) and Eq. © replacing the integration variable A with R; in this way 
we extract 

J &corrS v NdR 



0.009. 



(13) 



J S v pfdR 

The corresponding AS amin — £ CO rr — allows us to answer to the question 
ii); in fact, as it is shown in Figure 1101 the vanishing point of AS amin is, 
according to the two -spin Ising model, nearby R rs 2 A that corresponds to 
the equilibrium configuration of the Hydrogen molecule. 




Figure 10: AS amin for the H2 molecule as a function of nucleus-nucleus distance. 



4 Differences between the Configuration Inter- 
action approach and the two— spin Ising model 



The model proposed in Sec. [T]provides us with a measurement of entanglement: 
indeed, Eq. ([3]) describes the von Neumann entropy as a function of coupling 
constant A, for small A. By using Eq. ([7]), we can express A in terms of corre- 
lation energy and substituting it in Eq. ([3]) we can obtain the variation of S v n 
in terms of £ corr . 



{ £ corrLog( 2{£ £ ; o °; r [ +2) ) + {£ corr + 4)Lo g ( 2( £ — ; r + 4 2) ) } 

SvN = Te ToTT — a ' *- 14 ) 

(b corr + 2)LogA 
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In order to calculated the coefficient of proportionality among S v n and S corr 
we make an expansion of S v n for £ C orr — * (or equivalently for A — > 0) at the 
first order, obtaining a straight line characterized by an angular coefficient given 
by ms vN (£ corr ) = + Lo^)- Since this behavior is uncorrect to represent 

the logatithmic singularity of S v n in the origin, we make an expansion of Eq. 
(|14p , preserving the logarithmic deviation, and we obtain an expression of the 
form 

SvN — A£ corr -\- B£ corr L/Og(£corr) j (15) 

where A = 1/2 and B = -l/(4Log2). 




0.02 0.04 0.06 0.08 0.1 
Ecorr 



Figure 11: A comparison among the behavior of Eq. (|14p and its linear approx- 
imation and the logarithmic one, for the Ising model. 



In order to compare the behavior of S v n in Eq. (|14p . we have organized 
the numerical data, calculated with the method proposed in pQ, by making a 
correspondence between each value of £ CO rr and its respective value of S v n, 
obtaining the plot in Figure [12] 




0.01 0.02 0.03 0.04 0.05 

Figure 12: A correspondence of £ C orr and S v n calculated by the numerical 
procedure suggested by [1] 



Of particular significance is the fact that, in the range where S is monotoni- 
cally increasing, the correlation energy has its maximum, consequently S seems 
to be not a function. Moreover, it is important to note that £ corr begins to 
decrease for R > 1 A, region where the states become mixed, i. e. ,Trp ^ Trp 2 ; 
as depicted in Figure fT3l 
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1.4 
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Figure 13: The increasing of the degree of mixing in the two electron state: in 
black we depict the trace of p, in red the trace of p 2 . 



Probably, for this reason, the procedure adopted in [Tj seems to be not cor- 
rect: the density matrix, in fact, is calculated starting by the excitation coeffi- 
cient of a wave function obtained developping with the Configuration Interaction 
Single Double method a pure two electrons state. 

However, even if we consider only the first branch of the plot in Figure [T2l 
i.e. , the numerical values of S v n corresponding with increasing values of £ CO rr, 
and we fit the values around £ CO rr — ► with a F = A£ corr + B£ corr Log(£ CO rr) 
we draw out numerical values of the coefficient different from the ones used in 
Eq. (fT5|) . This result is shown in Figure [TJ] 




0.01 0.02 0.03 0.04 0.05 



Figure 14: A fit of S v n as a function of £ corr , around the origin, with a function 
of the form F = A£ corr + B£ corr Log(£ corr ) whose coefficients A and B assume 
the numerical values in Figure. 

In particular the arithmetic sign of the coefficient B in the two models are 
opposite and this implies the opposite concavity of the curve. 

This fact, clearly demonstrates a not satisfactory agreement between the 
Ising model and the one proposed in pQ. 
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5 Concluding remarks 



We have explored the role of entanglement in the model of two qubits describing 
the Hydrogen molecule JT]), considered as a bipartite system. In our discussion 
we have limited to the ferromagnetic case governed by the interaction coupling 
parameter J > 0. 

The concept of entanglement gives a physical meaning to the electron cor- 
relation energy in structures of interacting electrons. The entanglement can be 
measured by using the von Neumann entropy or, alternatively, the notion of 
Concurrence [7]. To compute the entanglement it is convenient to consider two 
Regions, say I and II, which provide two different reduced density matrices. 
The entropy turns out to be an increasing function of the coupling constant A 
in Region I, but the state under consideration is maximally entangled in Region 
II indipendently from the anisotropy parameter g. 

An interesting result is that for large coupling constants the entropy ap- 
proach 1, meaning that all levels are equiprobably visited by the considered 
spin. 

For weak interactions, at the boundary point At, = . 2 the von Neumann 

entropy admits a discontinuity, indicating a crossing of the lowest eigenvalues 
and, in a more general constext, a quantum phase transition [5]. 

In Sec. [5] a comparison between the entanglement and the correlation energy 
is performed. 

To quantifying the entanglement we resort to the reduced density matrix. 
The entanglement can also be measured by exploiting the concept of Concur- 
rence. 

The entanglement measure is always bounded, while the energy correlation, 
Zcorr — |£b| — 2 = \/4 + A 2 — 2, is a divergent function of A. This fact tells us 
that to look for simple relations valid on the whole A— axes has no sense. 

Thus, by limiting ourselves to weak couplings, we have minimized the mean 
square deviation given by Eq. (jHJ). This procedure leads to the value a m i n » 
—0.691217 for the minimizing parameter (see Eq. ©). 

Sec. [1] contains a comparison between the von Neumann entropy and the 
Concurrence. 

Such a comparison is illustrated in Figure[Tl for two spin system as a function 
of the coupling A for g = 1. 

Some important points are commented in Remark 1 and Remark 2 . 

In Figure 2] the deviation and its derivatives with respect to A are 

computed and a m i„ is evaluated for a ranging in the interval —1.29 < a < 
-0.091. 

In Figure [5] the minimized Concurrence deviation AC («) for the four eigen- 

states of the 2-spin model is shown. 

We point out the existence of a perfect symmetry among the Concurrence 
deviations for pairs of eigenstates of opposite eigenvalues. 

Formula (jTTj) . due to Heitler-London [SJ, is reported, where the position 
dependent exchange coupling J(r) is expressed in term of the length r of the 
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nucleus-nucleus separation in the Hydrogen molecule. 

To conclude, the magnetic field B has been monitored such that the equi- 
librium configuration more closest to the experimental one, r « 2.00, is the 
minimum occurring at r = 1.88 for B — | Ry and r — 2.25 for B = 0.375 Ry 
for the function AC a > . (r). 

We observe also that in the intermediate range of values, i. e., for 0.6 < A < 
1, the two functions SvN and the correlation energy are almost proportional 
within the 10%. 

However, when we organized the pairs of points (£ C orr, S v n) calculated by 
following the procedure described by [1], it is clear that the von Neumann en- 
tropy cannot be considered a function of correlation energy. The principle cause 
is that the function E corr presents a maximum in the region where S v n is mono- 
tonically increasing. 

The reversing behavior of correlation energy occurs in correspondence with 
an increase of the mixing degree of the two electrons state. The function £ CO rr 
in terms of the nucleus - nucleus distance R, increases till the state is pure, on 
the contrary, when Tr(p 2 ) becomes discordant from Tr(p), the function £ corr 
decreases. 

This fact suggests us that the numerical model based on the calculation of 
S v n starting by the excitation coefficients Ci, isn't completley correct because 
the density matrix is obtained as a product of two electron pure states. However, 
even if we consider only a branch of the plot in Figure [T2l the function obtained 
by the two spin Ising model, i. e., Eq. (| 14[> . is unsuitable for fitting these 
numerical data. 

On the basis of our results, essentially grounded on numerical considerations, 
in the near feature we would explore more complicated systems of molecules, 
such as for example the ethylene or other hydrocarbons, and compare these 
studies with the goals obtained for the Hydrogen molecule. 
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